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I I nucleate at very large distance scales, they can grow as detonations as well as deflagra- 
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1 Introduction 



Much of the interest in the cosmological first-order QCD phase transition has been 
due to the possibihty of affecting the big-bang nucleosynthesis [|[] . This would require 
creating significant inhomogeneities in the baryon number density, with a character- 
istic distance scale of the order of 10~Hh. The typical distance /„ between bubbles 
nucleated in the phase transition is related to the bubble surface tension, /„ oc a^^"^. 
The excitement has diminished as lattice calculations have yielded rather low estimates 
of a, indicating much smaller distance scales. 

The purpose of this paper is to point out that a low a does not necessarily imply 
short distances. The reason is that in addition to a, the distance scale depends on the 
latent heat L. As the phase transition appears to be at most only weakly first order, 
surface tension and latent heat should both be small compared to energy scales at the 
transition. Using the bag equation of state and classical nucleation theory, we find 
that even for small a, there is a large parameter space where the phase transition is 
preceded by considerable supercooling, and the distances between critical bubbles are 
large. Interestingly, we also find that with a more realistic equation of state than that 
of the bag model, the distances get even larger. It will be very enlightening to redo our 
calculations, when the whole equation of state of the QCD matter can be extracted 
from lattice calculations. 

In view of nucleosynthesis, not only the distances between critical bubbles but 
also the later stages of the phase transition are important. Since we do not know the 
microscopic physics operating at the phase transition front, we cannot give a definite 
solution to the problem. However, we can analyze what processes are allowed by hydro- 
dynamics. We find that the hadron bubbles could grow as detonations and that there 
is not necessarily a stage of slow growth in thermodynamical equilibrium during the 
phase transition. This leads to a rather unconventional picture of the phase transition. 
We also find that if this picture is correct, then it is unlikely that the QCD phase tran- 
sition has a noticeable effect on nucleosynthesis, in spite of the large distance scales. 
On the other hand, it is also possible that the hadron bubbles grow as deflagrations, the 
distance scales are equally large as above, and the QCD phase transition may have a 
noticeable effect on nucleosynthesis. Clearly, for conclusive results, a better knowledge 
of the parameter values must be acquired. 

The plan of the paper is the following. In Sec. |2JNucleationsection.2 



we review 



classical nucleation theory and calculate the nucleation temperature and the distances 
between critical bubbles in the QCD phase transition, using the bag equation of state. 
In Sec. |3The effect of the equation of statesection.3| we discuss the effects of a more 
realistic equation of state. Sec. [4Detonations and defiagrationssection.4| contains a 
hydrodynamical analysis of the possible growth mechanisms of the hadron bubbles, 
and in Sec. ^Baryon numbersection.5| we discuss the effect of the QCD phase transition 
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on nucleosynthesis. The conclusions are in Sec. |6 Conclusionssect ion . 6 



Since the hydrodynamics of bubble growth in this context has sometimes been 
analyzed using a plane- symmetric geometry {i.e., 1+1 dimensions), we stress that in 
this paper we use spherical geometry everywhere. 



2 Nucleation 

We review the classical thermal nucleation of thin-wall bubbles with the bag equation 
of state ^ 

Pg = QqT* - -, Ph = ahT^, (1) 

where L = 4(aq — ah)T^ is the latent heat. We use the notation T = T/Tc, where Tc 
is the critical temperature of the phase transition. The radius of a critical bubble at 
temperature T < 1 is 

2. 8. 1 



Ph-Pq L l-T^ 

where a is the bubble surface tension. The free energy expended to nucleate a critical 
bubble is 

N 4:71 9 167r cr^ , , 

Win) ^ -arl ^ -J^r (3) 

The nucleation rate is 

pit) = CT^-^W, (4) 

where 

Sit) ^ Wir,)/T = \ , (5) 

and C is a prefactor roughly of order unity. 

There are thus two essential parameters, L and a. The latent heat L gives the 
difference of and ah. Otherwise the values of and affect the nucleation only 
slightly through the expansion rate of the Universe. We are interested in the case of 
a small L. Therefore, and ah cannot be taken to have their ideal gas values — 
51.257r^/90 and 17.257r^/90, which correspond to a large latent heat L = 14.9T^ — 
but must lie closer to each other somewhere between these values. In reality they are 
functions of temperature ag(T), a/i(T), but until Sec. |3'l'he effect of the equation ol 



tatesection.3| , we take them to be constant. 

With constant and small L, the expansion (cooling) timescale is approximated 

by 

tT^ = const. = t^T^. (6) 

The Universe reaches the critical temperature at t = tc. The Universe then supercools 
to T < Tc and the nucleation rate begins to increase rapidly. A nucleated bubble 
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grows by detonation or deflagration, preventing further nucleation inside the volume 
(47r/3)t'^(t — t')^, where t' is the nucleation time. In the case of detonation the relevant 
velocity v is that of the phase boundary (detonation front), v = v^ct- In the case of 
deflagration, it is that of the shock driven ahead of the phase boundary, v = Vsh- In 
both cases, < f < 1, where = l/y/S is the speed of sound. The expansion rate 
of the Universe can be ignored during this short period of rapid bubble growth. The 
fraction of space affected is then 

ft 47r 
/ dt'p{t')—v^{t-t 
Jtc 3 



F(t) = 1 - exp 



(7) 



This grows rapidly from F ~ to F ^ 1 at the phase transition time t/, which we 
define by 

F{tj) ^ 1 - 1/e. (8) 

Bubble nucleation then ceases. 

Let us note that for very slow deflagrations, fdefi~0.1, the shock preceding the 
deflagration front is extremely weak. Between the shock and the deflagration front, 
the temperature increases continuously [^. However, even the temperature at the 
deflagration front, T^, may have been raised so little above the nucleation temperature, 
Tf, that the nucleation of new bubbles is not appreciably suppressed. Thus the velocity 
V should be the deflagration velocity Udefl, instead of the shock velocity v^h- As will be 
seen below, this leads to much shorter distance scales. 

Approximat ing 

S(t)^S{tf)-S'{tf){tf-t) (9) 
and extending the integral to — oo, Eq. ^ becomes 

S'itfY = 87Tv^CT^e-^^'f\ (10) 

The typical distance /„ = n~^^^ between centers of neighboring bubbles is obtained 
from their number density 

(11) 



The above approximations require 



\S" 



< 1 (12) 



and 



From Eqs. ^ and ^ we obtain 



\S\ 

|5'|(t/-g>i. (13) 



Sitf) = (14) 
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and 

Aff 



where 



S"{tf) = 'iii,(63ff + 3iff-l), (16) 



A^'-^^ (IT) 



3 L^T, 
and 

y^l-ff. (18) 

The nucleation time tf and temperature Tf are solved from Eq. 0. Taking the loga- 
rithm and using Eqs. |14| and ^ it becomes 

y' = = , (19) 

TfS{tf) Tf[SA-^ln\S'{tf%\y 
where 5*^ = InSTTv^CT^t"^ ~ 4 In T^tc — 170. This gives an equation for y, 

^1/2 



y 



;i - y) V8 [5^ - 4 In I A(l - y) 1/4(8 - 9i/)/(22/3) |1 ' 



(20) 



which can be solved iteratively. The phase transition temperature is then 

Tj = {l-yY/'T, (21) 

and the typical bubble distance 

= ''"l(l-,)v'(8-9,)«'" 



where tn = '^tc is the Hubble distance at t = tc. We plot these quantities in Fig. |1 Conclusionsflgur~ 
They depend on the parameters a and L only through the combination |1^. 

The approximations ^ and |T3] fail near / L'^Tc = 0.25 when the supercooling 
1 — Tf and the distance scale In are becoming very large. The nucleation rate, according 
to Eq. 1^, is beginning to grow more slowly and has a maximum at T = Tc/ ^/S. This 
would indeed lead to extremely deep supercooling and very large distances. Of course, 
our equation of state |l], tuned to be valid near T^, is then no longer applicable. 

The critical radius 

2a 8a , , 

Ph - Pq Ly 

must also not be too small (~ T~^) for the above thin- wall nucleation to apply. From 
Fig. |2Conclusionsfigure.2| we see that for small a, the critical radius is reasonably large, 
and almost independent of L. 
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The values of L and a indicated by pure glue lattice Monte Carlo simulations are 
of the order L ^ 2T^, a ~ 0.02T^^ |], ^. However, the uncertainty in these values is 
very large and, in particular, there is no lower limit. Thus, based on present knowledge, 
either one could be arbitrarily small. From Fig. |2Conclusionsfigure.2| we see that, e.g.^ 
length scales of In/vtu ^ 10~^-10~^ are possible. This is a distance scale large enough 
to affect nucleosynthesis. The corresponding critical radii Tc are large enough so that 
the thin-wall calculation should be valid. 



3 The effect of the equation of state 

Lattice Monte Carlo simulations imply that the energy density must have a very strong 
variation within a narrow temperature interval (< 10 MeV) in QCD with physical quark 
masses 0. Combined with the smallness of the latent heat this means that the realistic 
equation of state must differ even qualitatively from that of the bag model. This is 
illustrated in Fig. |3Conclusionsfigure.3| . As has been discussed above, the parameter 



values of the naive bag model can be corrected to reproduce a desired latent heat. 
However, this corrected bag model ( "bag model of text" ) cannot mimic at all the strong 
variation in the energy density. In this section we will analyze how the large derivative 
of the energy density curve, i.e., the large heat capacity, affects the nucleation. 

In the case of the nucleation calculation we may assume that the Universe is, at 
least locally, in the quark phase and very near thermal equilibrium. Furthermore, we 
can ignore the tiny quark and lepton chemical potentials. We denote the cosmic scale 
factor of the Robertson- Walker metric by i?, the energy density^ by e, and the entropy 
density by s. From the equation d(ei?^) = — pd(-R^) it follows that 

R 3(e + p) 3s(T) T ' ^ ' 

which can also be written as the equation d(si?^) = for entropy conservation. Here 
c(T) is the density of heat capacity (specific heat) of the quark phase in constant 
volume. It satisfies the relations c(T) = de/dT = T(ds/dT) = T(d^p/d^T), where p 



is the pressure. Equation |2J tells how the relation between expansion and cooling 
depends on the equation of state. Notice that in the case of the bag equation of state 
the factor c/3s equals unity, and the simple relation T oc 1/R is valid. 

The expansion rate of the Universe is determined from the Friedmann equation, 

^ - (2^) 
i?2 - 3M^, ' 



"'^In this section, all quantities are measured in the supercooling quark phase. 



5 



where Mp\ denotes the Planck mass. The vanishingly small curvature term was omitted. 



Eliminating dR/R from Eqs. ^ and gives the cooling rate of the Universe: 



3s(T) 87ie(T) T , 

If the energy density curve e(T) is very steep, the first factor is much smaller than 
unity. The Universe expands and the energy density decreases, but this causes only 
very slow cooling. In other words, the total energy density determines by Eq. ^ the 
expansion rate of the Universe, and this in turn determines by Eq. ^ the value of the 
derivative e(t); but then dT/dt = e(t)/c(T) is very small for large heat capacity. 

Since the Universe cools slowly, it has time to nucleate even at a large action, and 
therefore the true nucleation temperature is slightly higher than that of the bag model. 
To see this, notice first that compared with the time spent in the supercooled state, 
the period of time during which essentially all the nucleation takes place is short, and 



the cooling rate ^ is practically constant during this period. We define the parameter 

S = (27) 



where Tf is the temperature at the nucleation time tf. In Fig. |3Conclusionsfigure.3| , 6 
is roughly the ratio of slopes of the reahstic equation of state, and the bag equation of 
state. The nucleation temperature is determined from Eq. |T0|. However, now the pres- 
sure difference Ph(T) —pq(T) is unknown, and thus is also the nucleation action S(T). 
Still, we can find out what is the effect of 6: in the logarithm of Eq. |T0| the value of 
the constant S'a — 170 increases by 4 log 5 . As long as this increase is small compared 
with 170 the amount of supercooling, defined as 1 — Tf/Tc, decreases only little. 

Although a large value of 6 barely changes Tf, it has another significant effect. 
Physically the most important quantity related to nucleation is the average distance 
between nucleation centers immediately after the phase transition, defined between 



Eqs. |10| and |Tl|. Its dependence on 6 is seen from Eq. |Tl| to be 

In oc 6, (28) 

where the small correction coming from the change in Tf was left out. This result tells 
that a steep drop in the energy density e(T) increases significantly the typical distance 
between bubbles. This happens because when the cooling rate is lower it takes more 
time for the nucleation action S{T) to decrease by a certain amount. 

To get an order of magnitude estimate for 6, we assume that the energy density 
changes from its asymptotic value in the quark phase to its asymptotic value in the 
hadron phase within 10 MeV. This is an upper limit indicated by lattice calculations . 
Then we get 5 ~ 8. Thus the distance scales would be an order of magnitude larger 
than those shown in Fig. |2Conclusionsfigure.2|. If the true temperature interval for the 
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rapid change in energy density were smaller, say 1 MeV, then distance scales would be 
correspondingly larger, (5 ~ 75. 

Finally, we note that the validity conditions of the nucleation calculation, Eqs. |1^ 
and [13], do not depend in practice on 6 unless the specific heat of the quark phase, 
c(T), changes rapidly within the temperature interval Tc. . .Tf. 



4 Detonations and deflagrations 

A bubble of h phase surrounded by supercooled q phase has two modes of growth 
available, detonations and defiagrations [§, |lOl- In a defiagration bubble the fiuid 
inside the bubble is at rest, but the growing bubble is surrounded by a shock wave 
moving out ahead of the phase transition (defiagration) front. In a detonation bubble 
the phase transition (detonation) front advances into the fiuid which is at rest, but is 
followed by a rarefaction wave where the fiuid fiows outwards following the detonation 



front. Relativistic detonation bubbles have been discussed in Ref. and relativistic 
defiagration bubbles in Ref. 0. 

Consider fiuid fiow through the phase transition front in the rest frame of the front. 
The infiow is subsonic for a defiagration, Vi^ < Cs, but supersonic for a detonation, 
fin > Cs- These processes are further divided into weak and strong depending on the 
outfiow velocity. For weak processes the nature of the fiow velocity does not change, 
i.e., fout < Cs for weak defiagrations, and fout > for weak detonations. For strong 
processes it changes, i.e., Vont > Cs for strong defiagrations, and Vont < Cs for strong 
detonations. The case where tiout = Cs is called a Jouguet process. 

In the rest frame of the unaffected fiuid, all detonations and strong defiagrations 
move faster than sound, whereas weak defiagrations are subsonic. 

Because of restrictions on the fiuid fiow pattern from the bubble geometry, strong 



detonations are not possible for phase transition bubbles [0. (This is true also for 
1-dimensional "bubbles"). 

In classical combustion theory for chemical burning the internal structure of the 
combustion front rules out strong defiagrations and weak detonations P, |^. In par- 
ticular, the internal structure of a detonation front consists of a shock heating up the 
medium to initiate combustion, immediately followed by a defiagration. For a weak 
detonation, this defiagration would be a strong one. Thus the impossibility of strong 
defiagrations implies the impossibility of weak detonations. 

The internal structure of a phase transition front is different from a combustion 
front. Heating by a shock does not facilitate the phase transition, and the structure 
of a detonation front is not shock+defiagration. Therefore weak detonations are not 
ruled out |12|. Thus slowly growing bubbles {v < Cg) will be weak defiagrations and 
fast ones {v > Cg) will probably be weak detonations [T^. Strong defiagrations might 
also be possible in some cases. 
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We denote by Th and Tq the temperatures of the two phases at the phase transition 
front. For a detonation bubble, Tq will be the phase transition temperature derived in 
Sec. |2Nucleationsection.2| , Tq = Tf, but Th > Tf, and we can even have Th > T^. The 
rarefaction wave cools the h phase, so the final temperature will be below Th (and T^) . 
For a deflagration bubble, Th will be the final temperature of the h phase, whereas Tq 
is not the initial temperature Tf, but Tq > Tf (and could exceed Tc), as the q phase 
has been heated by the shock wave. 

The detonation and deflagration solutions are obtained from the hydrodynamical 
conditions of energy and momentum conservation. These processes must also satisfy 
the condition of non-negative entropy production [|], ^ These constraints do 

not fix the process uniquely for a given initial temperature Tf. Instead, we have a 
one-dimensional family of allowed solutions for each Tf, with different temperatures 
{Th,Tq) and different bubble growth velocities fdofi or v^et- This family may contain 
both deflagrations and detonations. Weak deflagrations are allowed for any Tf < 
but detonations and strong deflagrations require a minimum amount of supercooling, 
see Fig. [4 Conclusionshgure . 4 . 

Below we give exact results for the bag equation of state If the energy densities 
are scaled by the bag constant B = L/A = {aq — ah)T^, and the temperatures by T^, 
these results can be given in terms of a single parameter 



ah 



1 - 



> 1. (29) 



We identify a process by a point on the {eh/B,eq/B)-p\a.iae, where eh and are the 
energy densities of the two phases at the phase transition front. There are a number 
of special points in this plane (Fig. [4Conclusionsfigure.4D . 

Point C corresponds to Th = Tq = Tc. This is the limit of weak deflagrations as 
Tf — > Tc. It is at the fdcfi — ^ limit (the diagonal line through C). Point D is the 
(weak) detonation which requires the least supercooling. This point is at the Vdet ~^ 1 
limit (the diagonal line through D). The velocities change steeply near these diagonals. 
Points J and G are the Jouguet detonations and deflagrations, respectively, requiring 
the least supercooling. The coordinates {eh/B,eq/B) of these points are 

_ / 3 4r - 1 \ 

fl = fJr^,i±|Y (31) 



,r — 1 r — 1, 

J={sj + tj,sj-tj), (32) 
G=isG + tG,SG-tG), (33) 

where 

2t 

sj = -(1 + COS2/?), (34) 
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/ 2r 

tj = 2J -cos/? (35) 

V r — 1 



So = -(1 + 2 sin /3 + 2^3 cos /5 sin/?), (36) 

r — 1 



/ 2r 

tG = -\ -(COS/3 + VSsin/?), (37) 



and 



1 r + 1 

6 = - arctan \ / . (3J 

^ 3 Vr-l ^ 



Converting the energy densities to temperatures by 

rp4 _ ''^ ~ ^ rji4 rpA _ ''^ ~ ^ ~ rp4 fQn\ 

^h — o ^J-ci J-q—^ J-ci l-jyj 



we have the following results. The maximum temperature for which the weak detona- 
tions are allowed is 

Tf = Tg(D) = r-i/"Te. (40) 

At the detonation front we then have Tq = r'^^^Tc-, = r^/^Tc- The maximum 
temperature for which Jouguet detonations are allowed, is 

Tf = T,(J) 



r + 1 2 ^ 2^/r-l 

+ - cos 2/3 - - V2W cos P 



3r 3 3 



1/4 



Tc. (41) 



We always have Tg(J) < Tg(D) < T^. 

For deflagration bubbles Tg ^ Tf. To relate these two temperatures we have 
numerically integrated the flow equations for the region between the shock and the 
deflagration fronts 0. We denote the maximum temperature at which strong defla- 
grations are allowed by Tf{G). The temperatures Tg(D), Tg(J), and Tf{G) are plotted 
as a function of r in Fig. |5Conclusionsfigure . 5| . 

If the Universe supercools very much, the phase transition is not able to reheat it 
back to Tc. The limiting supercooling temperature is obtained from 

e,(r,) = e,(T,). (42) 

For the bag equation of state |I|, this gives 

T, = (— — ) T,= -] T,. (43) 

\ 3r / \r 3r J 

If Tf > Tj., the phase transition reheats the Universe to Tc before completing. The 
rapid detonation/deflagration stage is followed by a slower stage, where both phases 
coexist in thermal equilibrium at T^, and the phase transition proceeds only as the 
Universe expands [Q. If Tf < Tr, the detonation/deflagration takes the transition to 
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completion. If r > 4, i.e., B > eh(Tc), then reheating to Tc is possible from arbitrarily 
low temperatures. 

Comparing Eq. ^ with Eq. ^ we note that is always below Tg{D), although 
for small r, these are close to each other. Thus in those cases where the Universe has 
supercooled so much that it will not reheat back to T^, weak detonations are always 
allowed. The similar statement for Jouguet detonations becomes true for r > 1.644. 

In Sec. |2JN ucleationsection. 2| we related Tf to a and L. The parameter r depends 
on aq in addition to L. By making some assumption about a^, we can convert a and 
L to Tf and r, and classify points in the (L, a) parameter space according to which 
processes are allowed. This is done in Fig. |6 Conclusionsfigur e . 6 . 

As a concrete example, let us inspect the case a = 0.01. Now it is seen from 
Fig. 5Conclusionsfigure.6 that if L = 1, only weak deflagrations are allowed and the 
Universe does not reheat to the critical temperature. If L = 0.1, weak deflagrations 
are possible as well, and the Universe reheats back to T^. If L = 0.01 even Jouguet 
processes and strong deflagrations are allowed. 

To summarize the results of this section, for a given phase transition temperature 
Tf < Tc, there is a one-dimensional family of allowed bubble growth processes. This 
family will always include weak deflagrations. It may also include weak detonations 
(if Tf < Tg(D)), Jouguet detonations (if Tf < Tg(J)), Jouguet deflagrations (if Tf < 
Tf{G)), and strong deflagrations (if Tf < Tf{G)). Which of these allowed processes 
actually occurs, depends on the dissipative mechanisms internal to the front which 



determine the propagation speed of the transition p|, |T3|. Even though external 
conditions may allow detonations, the actual process could still be a slowly propagating 
weak deflagration. 



5 Baryon number 

Much of the interest in the QCD phase transition in cosmology stems from the possi- 
bility of leaving behind strong inhomogeneities in the baryon number, and maybe thus 
affecting big bang nucleosynthesis. The baryon number in the q phase is carried by 
massless quarks, but in the h phase it is carried by nucleons, with m ^ T^. The baryon 
number does not penetrate the phase boundary easily, and accumulates as a layer on 
the q side of the phase boundary. As more baryon number accumulates onto this layer, 
more will also leak through, but the net effect is that of dragging baryons towards the 
regions which remained longest in the q phase \W, ^ |r^, [T^. 



Even assuming we know the hydrodynamic details of the phase transition discussed 
above, it is difficult to estimate the shape and density contrast of the baryon number 
inhomogeneity. It depends on the rate of baryon transport within each phase and 
baryon penetration of the boundary, which are not known. 

To have a significant effect on nucleosynthesis, a number of conditions need to 
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be satisfied. 1) The distance scale should be large, > 1 m 10~Hh- 2) The 
high-density regions should contain most of the total baryon number. 3) The density 
contrast R, i.e., the ratio of baryon number density in the high-density region to that 
in the low-density region, must be large, i?finai > (?^p/?^n)nucicosynthesis ~ 7. 

We have argued in Sections |2Nucleationsection.2| and |3'l'he effect of the equation 
of statesection.3l that condition 1 could be satisfied. Because of the difficulty of getting 
baryon number through the phase boundary, condition 3 does not appear unreasonable. 
Condition 2 is perhaps the most difficult 

During the phase transition, baryon number has been collected onto a layer on 
the surface of the bubble, with some thickness d. The baryon density in this layer is 
-Riayer timcs larger than elsewhere. To have most of the baryon number in this layer, 
we must have roughly -Riayer > In/d. If the thickness of the layer is due to microscopic 
diffusion of baryon number away from the boundary (in the q phase), d will be very 
small. It has been argued that turbulent transport will be much more effective than 
microscopic diffusion, and lead to a much thicker layer, so that condition 2 might be 
satisfied |I9 |. 

In the usual picture of this phase transition, most of the growth of the h regions 
will happen in the equilibrium stage of the transition, very slowly. In the picture we 
have presented here, which leads to large distance scales /„ even with a small surface 
tension a, this stage does not usually exist. The phase transition is completed by 
detonation or deflagration. 

Especially, if the bubble growth process is a detonation, the phase transition is 
completed rapidly. The detonation bubbles grow with v^et > Cg. In those regions 
where the detonation bubbles colhde, the transition is then already completed, and the 
turbulence caused by the collision will operate only in the h phase. Where there is any 
q phase left, the detonation front will move on unaffected, and once they have covered 
the space between them, the transition is over. The baryon number accumulating 
on the phase boundary cannot escape from the supersonically moving front. Thus the 
layer should get no thicker than a few fm. This appears to make condition 2 impossible 
to satisfy. 

If the bubble growth process is a deflagration, the phase boundary will be slower 
and turbulent baryon transport may be effective in making the layers thicker. As 
mentioned in Sec. pNucleationsection.2| , the distances between nucleated bubbles do 
not depend on whether the bubbles grow as deflagrations or detonations, unless the 
deflagrations are exceedingly slow. Therefore, slowly growing deflagration bubbles, 
with efficient turbulent baryon transport, seem to be able to affect nucleosynthesis. 

If the hadron bubbles grow as deflagrations, but with an exceedingly low velocity, 
there may be another mechanism, in addition to turbulence, for enhancing the baryon 
number. Neutrinos, which do not carry any baryon number, may carry a considerable 
part of the energy flux. The reason is that the hydrodynamical flux wy'^v, measured 
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in the rest frame of the phase transition front, vanishes as v 0. This curve is given 
as the solution of the equation Pq(Tq) = ph{Th). The solution does not agree with 
the curve Tg — Ty^^ so the neutrino flux 7r^/120)(T|^ — T^), where g^n denotes the 
effective number of active neutrino degrees of freedom, remains non-zero. But as noted 
above, for these very slow deflagrations the distance scales are also small. However, 
according to Sec. |3'l'he effect of the equation of statesection.3| an extremely large heat 
capacity would increase the distance scale with several orders of magnitude. There is 
also another escape: suppose that the initial growth mechanism is a deflagration which 
is not very slow, so that distance scales are large, but that then the Universe reheats 
to Tc, where a stage of very slow growth near thermodynamical equilibrium follows. 
Here, neutrino transport may be effective [^, in addition to turbulence. If this stage 
lasts long enough, the baryon number remaining at the high-density regions could be 
enhanced. From Figs. |2Conclusionsfigure.2| and |6 Conclusionsfigur e . 6| one can see that 
this scenario is possible in a very small but non-vanishing region of the parameter 
space, near e.g. the point L ^ 2T^, a ~ 0.3T^. 

Finally, let us note that so far we have only studied which hydrodynamical pro- 
cesses are in principle possible, without being able to fix a definite growth velocity. In 
Ref. [|T^ a model is presented which tries to fix this velocity, by introducing a phe- 
nomenological dissipative constant F. With a dimensional estimate F ^ 1 T~^, and the 
parameters L ^ O.IT^, a ~ O.lTj^, the bubbles grow as deflagrations, with a velocity 
Vdefi = 0.1. The neutrino flux is vanishingly small, but turbulence might be effective. 



6 Conclusions 

Parametrizing the QCD phase transition with the latent heat L and the surface ten- 
sion a, we have studied bubble nucleation and growth. The effects of a more general 
parametrization have been estimated. We have investigated the possibility that the 
inhomogeneities generated at the QCD phase transition significantly affect nucleosyn- 
thesis. 

We find that parameter values in the range of, e.g., a ~ 0.01-O.lT^^, L ^ 0.01- 
O.IT^"*, lead to a transition with large supercooling, relatively large critical bubbles, 
and a distance between bubbles of ~ 10~^-10~^tff. These bubbles may grow as 
detonations as well as deflagrations, and the Universe does not reheat to Tc. In spite of 
the large distance scale of the inhomogeneity, phase we found that for detonations the 
accumulation of baryon number in the high density regions is too ineffective to make a 
noticeable effect on nucleosynthesis likely. Even in the new region of parameter space 
studied, only deflagrations with efficient turbulent baryon number transport, or maybe 
with large neutrino flux, seem able to affect nucleosynthesis. 
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FIGURE CAPTIONS 



Fig. 1: The supercooling 1 — 7/ and the bubble distance scale In- The distance scale is 
given as l„/vtH, where tn is the Hubble distance ("horizon") and v is the detonation 
or shock velocity, l/-\/3 < v < 1. The dependence on the surface tension a and 

the latent heat L is through the combination a^/L^T^. We also show the nucleation 
action S{tf). The other thin lines show the two quantities S" {t f) / S' {t (short-dashed 
line) and [|5"(t/)|(t/ — tc)] ^ (dotted line), whose smallness we have assumed. Our 
approximations are seen to break down at cr^/L^Tc ^ 0.25. 

Fig. 2: Contours of the critical radius Vc and the bubble separation Z„ on the (L, a) 
parameter plane. The sohd line corresponds to a^/LF'Tc — 0.25. This figure is for the 
bag equation of state, equation of As discussed in the text, use of a more realistic 
equation of state could increase the distances In by an order of magnitude or more. 

Fig. 3: Schematic representation of the energy density versus for three different 
equations of state. Thinner parts of the curves denote the metastable branches. For 
clarity the magnitude of L has been exaggerated in the figure. 
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Fig. 4: Detonations and deflagrations: This plot shows how the different processes lie 
in the {eh, eq)-plane. The entropy condition restricts the allowed processes below the 
AS — curve. Point C corresponds to Tq — Th — Tc- For a given Tf there is a 1- 
dimensional family of solutions, denoted by the dashed hne. The detonation branch of 
this family is a horizontal line, the deflagration branch a steep curve. For any Tf < 
it always passes to the left of point C, indicating that weak deflagrations are allowed. 
If it passes to the left of G, strong deflagrations are allowed. If the detonation branch 
passes below D (J), then weak (Jouguet) detonations are allowed. This figure is for 
r = 1.01. 

Fig. 5: Some special temperatures for the bag equation of state as a function of r. 
Tq{D) is the maximum temperature for which (weak) detonations are allowed. Tq(J) 
is the maximum temperature for which Jouguet detonations arc allowed. Tf{G) is the 
maximum temperature for which Jouguet deflagrations are allowed. For T < Ty(G), 
strong deflagrations arc allowed. is the lowest temperature for which the latent heat 
is sufficient to reheat the universe back to T^. 

Fig. 6: Regions on the (L, a) parameter space, where different processes are allowed. 
This figure is for a bag model with = 34.257^V90 + L/87;^ = 34.257rV90-L/87;l 
The solid lines divide the graph in three regions depending on what kind of detonation 
bubbles are allowed by hydrodynamic considerations. Weak deflagrations are always 
allowed. Strong deflagrations are allowed above the long-dashed line. The universe 
reheats to if wc arc below the short-dashed line. Thus detonations are always 
allowed in those cases where the universe does not reheat to Tc. 
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